The Effective Temperature in Elasto-Plasticity of Amorphous Solids 
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An effective temperature T e g which differs from the bath temperature is believed to play an es- 
sential role in the theory of elasto-plasticity of amorphous solids. The definition of a measurable T e g 
in the literature on sheared solids suffers however from being connected to a fluctuation-dissipation 
theorem which is correct only in equilibrium. Here we introduce a natural definition of T c s based 
on measurable structural features without recourse to any questionable assumption. The value of 
T c ff is connected, using theory and scaling concepts, to the flow stress and the mean energy that 
characterize the elasto- plastic flow. 



Introduction: Amorphous solids form when super- 
cooled liquids are further cooled below the glass transi- 
tion. While indistinguishable in their microscopic disor- 
der from fluids, amorphous solids exhibit, in contradis- 
tinction from fluids, a yields stress below which they re- 
spond clastically to external strains; fluids flow under any 
external strain. Once the amorphous solid is subject to 
large enough strains such that the response of the internal 
stress exceeds the yield stress, it can flow plastically in a 
manner that depends on the temperature, the shear rate, 
the density etc. While there have been many attempts to 
present phenomenological equations to describe the rhe- 
ology and the constitutive relations of such elasto-plastic 
flows [H, 0, H, H H, H, 0], to this date none of these at- 
tempts has gained universal acceptance. In fact, there is 
no complete agreement even on the field variables, or 'or- 
der parameters' which are necessary to close a complete 
set of equations. 

Among the more interesting ideas for order parame- 
ters stands the proposition that such elasto-plastic flows 
exhibit two different temperatures, the regular tempera- 
ture T that relates to the mean velocity of the particles 
forming the amorphous solid (and, naturally, to the heat 
bath to which the system is coupled), and an 'effective 
temperature' T e g that has to do with some 'noise' § or 
'compositional properties' of the material The con- 
cept of effective temperature was originally formulated in 
Ref. [10I | in order to describe the macroscopic properties 
of granular matter. Although these systems are ather- 
mal in the sense that they cannot evolve under normal 
temperature conditions, they may be driven by an exter- 
nal force which establishes an equilibrium like behavior 
described by an effective temperature related to the in- 
tensity of the "tapping". These ideas have since been 
applied in other kinds of mechanically driven systems 
such as disordered elastic structures 
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12]. For sheared 
amorphous solids the definition of this effective temper- 
ature was never entirely convincing, being connected ei- 
ther to the density of ill-defined 'shear-transformation 
zones' 0, Jsl or related to some 'fluctuation-dissipation 
theorem' [13L HH, EH that should be exact only in equi- 
librium [1 61 ] . The aim of this Letter is to announce a sim- 
ulational discovery of very sharp and direct meaning to 
an effective temperature in a number of simple computer 



models of elasto-plasticity. This effective temperature 
has an obvious connection to the compositional disorder 
in the material. Moreover, it naturally identifies with 
the regular equilibrium temperature in the super cooled 
liquids. In this way it allows a smooth description that 
unites the super-cooled regime with the amorphous-solid 
regime, something that is certainly lacking in many phe- 
nomenological descriptions. To define T e g we need first 
to recall some recent advances in describing the super- 
cooled equilibrium regime. 

Up-scaling in the super-cooled regime: In a se- 
ries of recent papers (cf. [llj and in particular [l8j]) it 
was proposed that the scenario of the glass transition, 
including the astonishingly rapid slowing down of the dy- 
namics in a short range of temperatures, is usefully en- 
coded by the temperature dependence of the concentra- 
tions of a finite set of quasi-species which can be indexed 
by 1, 2, • • • ,n. The precise nature of these quasi-species 
may change from model to model, but they are always 
formed by particles and their nearest neighbors. The 
main advantage of these quasi-species is that they obey 
a discrete statistical mechanics, in the sense that their 
temperature-dependent concentrations (Ci)(T), are de- 
termined by a set of degeneracies gi and enthalpies Hi 
such that 



<C;)(T) = 



9ie 



-Hi/k B T 



(1) 



Obviously, if such a simple description is available, we 
can predict which quasi-species will be there when the 
temperature is high, and which will remain when the 
temperature decreases: only those with lowest free en- 
ergy remain at low temperatures. Indeed, simulations 
show how the concentrations of some quasi-species de- 
crease, some increase, and some start increasing and then 
decrease when temperature is lowered, according to their 
degeneracy and enthalpy. Fluidity (or short relaxation 
times) is therefore correlated with high concentrations 
of quasi-species whose free energy is high, and solidity 
(or long relaxation times) is correlated with high concen- 
trations of quasi-species whose free energy is low. This 
qualitative observation was made quantitative by noting 
the concentrations of those quasi-species that tend to dis- 
appear when the temperature is lowered, and summing 
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these concentration to what was called the 'liquid-like' 
concentration (Cf)(T). The inverse of this concentra- 
tion provides a length-scale (the typical distance between 
'fluid' quasi-species): 

£(T) = [{Ct){T)]- l ' d , e^oowhenT^O, (2) 

where d is the space dimension. It was amply demon- 
strated on a large variety of models that the relaxation 
time T a (T) measured using correlation functions in the 
super-cooled regime is determined by this diverging scale 
according to 



(3) 



where \i is a typical free energy per particle and To a 
microscopic (cage) time. In contradistinction with the 
Vogel-Fulcher and Adam-Gibbs fits, Eqs. ^ and © 
imply that there is no singularity associated with the 
glass transition at any temperature other than T — 0, as 
was explained in fl9j | . The discovery that we announce 
here is that the crucial statistical mechanical relation Eq. 
(TTJ) can remain correct and very useful, with T replaced 
by T c g , also in the elasto-plastic regime of the amorphous 
solids that form at ultra-low bath temperatures. In other 
words, T e g exists and it determines the compositional 
disorder of the material. 

Two Models: We present the new findings using two 
different models of glass formation in two-dimensions, the 
first being the Shintani-Tanaka model [2(|, and the sec- 
ond the so-called 'hump model' which was inspired by 
[2~fl ] and analyzed in Ref. [Hl |. The Shintani-Tanaka 
model has N identical particles of mass to; each of the 
particles carries a unit vector itj that can rotate on 
the unit circle. The particles interact via the potential 
U(r tJ ,e u 6j) = U(n j )+AU(r lJ: e i , 6j). Here Ufa) is the 
standard isotropic Lennard- Jones 12-6 potential, whereas 
the anisotropic part AU(rij,8i,9j) is chosen such as to 
favor local organization of the unit vectors in a five-fold 
symmetry, to frustrate crystallization. For full details of 
this model the reader is referred to Refs. [2(il. I22I [23j ; here 
it suffices to know that with the parameters chosen in 
Ref. [20( the model crystallizes upon cooling for A < 0.6 
whereas for larger values of A the model exhibits all the 
standard features of the glass transition, including a spec- 
tacular slowing down of the decay of the correlation func- 
tions of the unit vectors Cfl(i) = (1/^V) W ' M i(0)) 
which is very well described by Eq. ([3]). 

The 'hump model' again employs N identical parti- 
cles interacting via a potential that is constructed as 
a piecewise function consisting of the repulsive part of 
a standard 12-6 Lennard-Jones potential connected at 
ro = 2 1 / 6 <7 to a pol yno mial interaction P(x) = a,iX l . 
The di's are tuned [18| so that P(x) displays a peak at 
r = TVimp and also such that there is a smooth conti- 
nuity (up to second derivatives) with the Lennard-Jones 
interaction at U(ro) — eho as well as with the cut-off in- 
teraction range f/(r*) — 0. The interaction potential for 
the hump model is shown in Fig. [T] Note that the two 




FIG. 1: The pair-wise potential for the hump model. In the 
inset we show a snapshot of the position of the point-particles 
(the circles represent the finite range of interaction, cf. 




FIG. 2: Upper panel: trajectories of T c s for the hump model 
as a function time for a fixed strain rate 7 = 10 -4 and 
N = 6400, settling at the steady state either from above or 
from below. Right lower panels; representative stress-strain 
curves for the seme model for different bath temperatures. 
Lower left panel: The measured steady-state concentration of 
the four relevant quasi-species of the hump model (full sym- 
bols), compared with predictions of Eq. ([l]) but with T a s 
(empty symbols). Circles, triangles, squares and pentagons 
correspond to (d) with i = 3, 4, 5 and 6 respectively. 



typical distances that are defined by this potential, i.e. 
the distance at the minimum r nim and the cutoff scale 
r*, appear explicitly in the amorphous arrangement of 
the particles in the supercooled liquid, as shown in the 
inset in Fig. [TJ The model has two crystalline ground 
states, one at high pressure with a hexagonal lattice and 
a lattice constant of the order of r m i n . At low pressure 
the ground state is a more open structure in which the 
distance r* appears periodically. At intermediate pres- 
sures the system fails to crystallize and forms a glass 
upon cooling [l8j . 

Simulations of the elasto-plastic regime: The 
molecular dynamics protocol we implemented for both 
models is the same. First, we carefully equilibrated a 
large number of independent configurations with TV par- 
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FIG. 3: Symbols: T c s a s function of Tb for the hump model at 
three values of the strain rate. Lines: prediction of T a ff using 
Eq. l(5). Inset: T e ff a s function of for the Shintani-Tanaka 
model at the two higher values of the strain rate. Straight 
lines indicate where T e s — Tt- 
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FIG. 4: Energy per particle in the elasto-plastic steady state 
as a function of the bath temperature. Symbols: simulation 
results. Lines through the symbols: theoretical prediction 
using Eq. Q. Black lines: equilibrium energy per particle. 



tides (N varied in these simulations between 1024 and 
6400) in the NVT ensemble using the Berendsen thermo- 
stat over a wide range of temperatures. These samples 
were then used to determine the enthalpies and degen- 
eracies for the two models as described in detail in 
These measured parameters predict accurately the con- 
centrations of quasi species at any given temperature as 
well as the r a relaxation time. After this, we turned our 
equilibrium super-cooled liquids into amorphous solids 
by minimizing their potential energy (conjugate gradient 
algorithm), allowing us to sample a representative set 
of meta-stable minima. This procedure can be thought 
of as quenching a liquid infinitely fast into a disordered 
solid whose temperature is formally T — 0. At this 
point we bring the particle velocities, using a short NVT 
run, to a value consistent with a desired bath tempera- 
ture Tfc. Then we force the system at a constant strain 
rate 7 using the SLLOD algorithm combined with Lees- 
Edwards boundary conditions. The stress-strain curves 
obtained for different bath temperatures are shown in 
the upper panel of Fig. for the hump model. After 
the usual elastic response, an irreversible plastic flow be- 
gins, eventually generating a time independent steady 
plastic flow in which all the thermodynamic quantities 
(such as flow-stress or energy) reached a value indepen- 
dent of the initial quenched configurations. The major 
finding is that throughout the evolution the concentra- 
tions of quasi-species obey Eq. (fTJ) with the equilibrium 
measured values of gi and Hi, but with T\> replaced by 
T c g which is a function of T), and 7. It is remarkable 
that there is an equally well defined effective tempera- 
ture not only in the steady state but also in the transient 
regime. The steady-state values of T e g as a function of T 
are shown for various values of 7 in Fig. [31 Note that at 
high temperatures T e ff — > T\, whereas T e g increases when 
Th — > 0, increasing the fluidity of the system. 

The relevance and importance of the steady-state ef- 
fective temperature can be demonstrated by relating it 
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FIG. 5: Flow stress as a function of Tt for the indicated values 
of the strain rate 7 for the hump model, and in the inset for 
the Shintani-Tanaka model for the two higher values of the 
strain-rate. 



to the mean energy and the flow stress, the latter being 
the mean stress in the elasto-plastic steady state. The 
mean energy per particle U/N was directly measured in 
the steady state of the hump model, and compared to 
the theoretical prediction 



U 
N 



0.835 



(4) 



where the first two terms were taken from the equilibrium 
theory for the hump model in Ref. [18j but with T e g re- 
placing Tb in determining the concentrations of the quasi- 
species; the last terms is the elastic energy per particles 
stored in the steady state. The almost perfect agreement 
between the theoretical expectation based on T c g and the 
direct measurement is demonstrated in Fig. [4] 

The dependence of the flow-stress on for various 
values of 7 is shown in Fig. We show now that we can 
predict the values of the flow stress cr 00 {T} > , A f) given the 
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FIG. 6: Double logarithmic test of the scaling function Eq. 
([5]); all the data from Figs. [3] and [5] are re-plotted here for 
both models (with the Shintani-Tanaka model in the inset), 
to demonstrate the excellent data collapse. The continuous 
line is the function f(x) — 4.93(x — 1) 



data for T e ff(Tb,7) or vice versa, predict T e g(T;,, 7) from 
the knowledge of Coo (Tt, 7). To this aim we invoke scaling 
concepts, and propose a scaling form for (Too (Tf,,T e s) such 
that the 7 dependence is carried here by T e ff: 

croo(Tb,T e s) = — k B T b f ( ) , (5) 

where p, m are the density and molecular weight of the 
particles. The function f(x) is a dimensionless scaling 
function that must obey /(l) = to agree with the ob- 



served fact that at higher temperatures where T e g = T b 
the flow stress approaches zero. For T — * or x — ► 00 
we observed that Goo becomes proportional to T e g re- 
quiring f(x) — > Cx for x — > co. The simplest function 
that obeys these limits is f(x) — C(x — 1). A test of 
the predicted data collapse is shown in Fig. [6] where the 
continuous line is the function 4.93(x — 1). Having the 
scaling function at hand we can predict the data, say of 
T c s(T b , 7) from the knowledge of tToo (T b , 7), or vice versa. 
Using for example the data for for the hump model 
in Fig. [5] and the scaling function f(x) = 4.93(x — 1) 
we solve for T e ff(X{,,7). The results are demonstrated in 
Fig. [3] with the curved lines going through the data. We 
conclude that the procedure is in satisfactory agreement 
with the data, demonstrating the importance of the con- 
cept of effective temperature. We mention in passing that 
the temperature T* where T e ff separates from T b can be 
easily predicted by equating the relaxation rate t q (£(T)) 
with 7" 1 : T* ~ ln~ 1 To7. At temperatures higher than 
T* the natural relaxation time r a is the shorter of the 
two, whereas at temperatures lower than the minimum 
the shear rate determines the rate of relaxation. 

Much remains to be done. For example in the case of 
the standard model of binary mixtures we did not find 
a satisfactory up-scaling that remains the same in and 
out of equilibrium. This and other riddles concerning 
the present approach will be dealt with in a future publi- 
cation. Nevertheless we believe that the present findings 
will provide grounds for new exciting future research. 

This work has been supported in part by the Israel 
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